Laser-like vibrational instability in rectifying molecular conductors 
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We study the damping of molecular vibrations due to electron-hole pair excitations in donor- 
acceptor (D- A) type molecular rectifiers. At finite voltage additional non-equilibrium electron- hole 
pair excitations involving both electrodes become possible, and contribute to the stimulated emission 
and absorption of phonons. We point out a generic mechanism for D-A molecules, where the 
stimulated emission can dominate beyond a certain voltage due to inverted position of the D and A 
quantum resonances. This leads to current-driven amplification (negative damping) of the phonons 
similar to laser-action. We investigate the effect in realistic molecular rectifier structures using first 
principles calculations. 

PACS numbers: 73.63.-b, 85.65.+h, 72.10.Di, 73.40. Gk 



In a seminal paper from 1974, Aviram and 
Ratner proposed an electronic rectifier based on 
a single organic molecule [1. Akin to the p-n 
junction in solid-state electronics their design in- 
volved a donor and an acceptor group bridged 
by a tunnel barrier (D-A). Ever since, the interest 
in molecule-based electronic operations, such as 
rectification [2]-[4], has increased |5]. Recently, focus 
has been on the interaction between the electronic 
current and the atomic dynamics in molecular 
conductors [6-8 . It is well known that a fraction 
of the electronic current exchange energy with the 
molecular vibrations leading to Joule heating [9- 
[I2J. However, instabilities such as conductance 
switching or break-down occur in experiments [H 
[T3UT6] in the high voltage regime, which still poses 
open questions to theory. Recent theoretical work 
has indicated that current-induced forces can con- 
stitute an important channel of energy exchange 
between electrons and ions different from Joule 
heating, and lead to destabilization or runaway be- 
havior of molecular contacts fTT} flS]. 

In this paper we point out a generic mecha- 
nism leading to instabilities in D-A molecular rec- 
tifiers. We show that the electron- hole pair ex- 
citation by the atomic vibrations, known as elec- 
tronic friction[19^, I20I, may not necessarily damp 
out the vibrational energy in biased D-A molec- 
ular rectifiers. Instead, we can get negative fric- 
tion or phonon amplification beyond a certain volt- 
age. This happens when the stimulated emission of 
phonon dominates over absorption processes due 
to a population inversion, similar to what happens 
in a laser. 

We first outline the theoretical basis, then illus- 
trate the effect by a simple two-level model, before 



presenting our first-principles calculations on real- 
istic systems relevant for experiments. 

Theory. — Within the harmonic approximation, 
ignoring mode-mode coupling, we can calculate 
the current-induced excitation of a given vibra- 
tional mode (phonon) using the rate equation ap- 
proach for the phonon population, N ^ 

N = B{N^l)-AN. (1) 

Here all phonon emission(absorption) processes 
are described by the rate, B{A)^ and we have the 
steady-state occupation. 



In equilibrium Nqs = nB{hCl), where is the 
Bose occupation at the phonon energy (fi(7), corre- 
sponding to ^ = e^^/^^^S > B. However, out of 
equilibrium, when an electronic current is present, 
we may obtain an occupation Nss ^ tib- Most 
importantly, we may obtain a 'laser'-type insta- 
bility, namely if B approaches A. Using Fermi's 
golden rule, we can calculate the rates resulting 
from the coupling to the electronic subsystem in 
the presence of current, 

B =^^\{f\M\i)\'F,{l-Ff)5{Aef.+fin) . (3) 
" ij 

Here M is the electron-phonon coupling between 
initial(i) and final(/) electronic states with occu- 
pations F, Asfi = Sf — Ei. A is given by the 
corresponding expression with Q ^ —Q. 

For an applied voltage, we employ the usual 
non-equilibrium setup with left(L) and right (i?) 
electrodes at different chemical potential (/iL — 



Mi? = eF), corresponding to Fermi distributions 
n^, a G {L^R}. We consider low temperatures, 
where ksT <C eV, fil]. In this limit B is only non- 
zero if eV > h£t^ and it can be written as, 

B= ALR{e,e-hn)de, (4) 

J iJ.R-\-hQ 

introducing the electrode-resolved, coupling- 
weighted DOS for the electron-hole pairs, 

Aa^{e,e') = ^TT[MA^{e)MA^{e')] , (5) 

given by the spectral densities, Aq,, of scattering 
states (DOS) originating in electrode a. In the 
same approximation the corresponding expression 
for A becomes, 

+ / ALL{e,£^hn)d£. (6) 
J ^L—hVt 

We shall now assume that we are dealing with a 
molecule with a donor level with energy £d which 
is most strongly coupled to the left lead, while an 
acceptor level with energy £a is primarily coupled 
to the right lead. In this case Alr{£^£') will have 
a resonance as a function of £ close io £d and a 
resonance as a function of £' around £A' The in- 
stability occurs when B > A. Inspecting Q, we 
see how B becomes large \i £d is in the energy 
window [jiR -\- HQ; fiL], ctnd £d — ^ £A' In this 
case A will be far from its optimum values since 
the integration intervals in (|6| are far from the 
resonances of Alr^ Arr^ All. The instability can 
also be described as the electron-hole pair damp- 
ing rate[19n21j, 

7^^ = ^ - S (7) 

going negative. In principle other damping mech- 
anisms such are coupling to bulk phonons, could 
be added to 7®^. 

Simple two-level model. — In Fig. [l^, we con- 
sider the simplest model with a donor- level (^i:) = 
—£q) which couple with the left electrode(ri:)), 
and an acceptor (s/^ = £q) coupling to right elec- 
trode(r^). If the DA hopping matrix element, t, is 
small, these levels will follow the nearby electrode 
chemical potential [22 . We introduce an electron- 
phonon coupling corresponding to modulation of 
the DA hopping 'distance', m cx dt/dx. 

The processes involved at T = 0, are shown 
in Fig. ^p-d. In equilibrium (/i/, = jiR = /io) 
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FIG. 1: (Color online) (a) Simple two-level model for 
the donor- acceptor (DA) system. (b-d)D- and A-DOS 
and filling(grey) at T = 0. The D(A)-levels follow the 
chemical potential, /j.l{ij^r) of the closest L(R) elec- 
trode, (b) For zero voltage eV = /hl — I^J^r = only 
absorption processes around 11 l — fJ^R = Mo are pos- 
sible. Similar excitations from the right electrode are 
not shown, (c) Absorption dominates for a voltage 
where sa — £d ^ and 7^^ is maximal, (d) Popula- 
tion inversion for larger voltage when, en — £a ^ hQ: 
Emission dominates and 7^^ < 0. (e) Contours of 
J^^(V) for varying t (dark is negative), £0 = 0.35eV, 
Ti = 0.3eV, Fs = O.leV, m = 0.005eV, HQ = 20meV, 
/io = -0.2eV. (f) 7^^(y) for t = 0.05eV showing the 
generic behavior (maximum followed by negative val- 
ues). 

electron-hole pair generation both intra- and inter- 
electrode lead to damping(absorption), see Fig.[T|3. 
The high D(A)-DOS around £d{£a) wih domi- 
nate processes when both DOS peaks are within 
the 'voltage window', [/j^r; /j^l]^ and have a differ- 
ent filling. As the bias is increased the £d{£a) 
moves up (down). A maximum damping is ex- 
pected when the filled D-states are separated from 
the empty A-states hy £a — £d ^ making ab- 
sorption (A) a dominating process, see Fig. 
More interestingly, as the bias is further increased 
so £d—£a ^ hftwe have a situation of 'population 
inversion': Electronic transitions from (filled) 
to 6:^ (empty) now makes phonon emission dom- 
inate over absorption, see Fig. [T]i. This corre- 
sponds to the Alr{£,£ — hQ)-teTm and thus B 
to be dominating. For fixed m and sufficiently 
small t we reach a situation where the friction be- 
comes negative when £a — £d ^ Fig. [l^. For 
big D — A coupling the effect disappears as the 
coupling makes bonding/anti-bonding DA-states 
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more relevant. For a particular choice of t the 
generic behavior of a maximum followed by nega- 
tive values of 7®^ is clearly seen, Fig. [T]F. 
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FIG. 2: (Color online) (a) Molecular STM junc- 
tion: S atom adsorbed on the left adatom("tip"), 
the right ("sample") adatom has adsorbed an NH2 
molecule, (b) Positions of the D {Spx^Spy of S) and 
A levels (2pz of N) as a function of bias. The dashed 
lines indicate /HL^fJ^n- Inset: Left/right DOS for IV 
bias, (c) Electron-hole damping(7^^) as a function of 
voltage bias for unstable phonon modes (inset) showing 
7"^^ < at high bias. 

We conclude that the phonon amplification in- 
stability, 7®^ < 0, is similar to a two- level atomic 
laser. In the presence of bias, the D state, filled by 
electrons, is located at an energy above the empty 
A state. In this case a phonon gets amplified by 
the stimulated emission accompanying electronic 
transitions from to sa- This is analogous to 
the population inversion and stimulated emission 
of photons in lasing. 

First-principles calculations. — We have per- 
formed NEGF-DFT based transport calculations 
[23l [24] on two types of experiments, where the 
instability appear. In the first example, we note 
that the tip in an STM can be functionalized by 
adatoms[25] or molecules [26l [27] enabling control 
of both voltage and hopping (~ t). This offers 
an ideal experimental test of the instability effect. 
In Fig. [2^ we consider a junction formed by two 
gold adatoms with S and NH2 adsorbed, repre- 
senting sample and functionalized tip. The Spx^y 
orbit als of S, located just below the jul, serves as 
D and follows jj^l- The 2pz orbital of N follows 
jiiR, and serves as A, see Fig. [2|b). In Fig. [2|c) 



we show the j^^{V) for the unstable modes with 
negative damping(insert). The behavior is simi- 
lar to the two-level model, and involve transverse 
adsorbate motion as selected by the symmetry 
of the D and A. In the next example the D-A 




^-O.T^^^^^-0.2 -0.1 0.1 0.2 0.3 



Energy (eV) 



FIG. 3: (Color online) (a)- (b) The A(LUMO) and 
D(HOMO) states at 0.3V (/xl = 0.15V,/iH = -0.15V 
shown as vertical lines) . (c) The position of D and A 
levels(DOS peaks) as a function of bias, (d) The Al, 
Ar at 0.3V. The dashed lines include only contribu- 
tion from the single D or A orbitals, while the solid 
lines include all contributions. 

behavior is tied to the surface anchoring groups 
of the 4-dimethylamino-4'-nitrobiphenyl molecule 
bridging Au(lll) electrodes [28l We iden- 

tify the D, A states from the molecular projected 
Hamiltonian[22] . The A-state(Fig. [3^) penetrat- 
ing to the left has a 7r-character, while the D- 
state with s-character(Fig. ^jp) exists near the 
left electrode. Again Al{Aji) mainly involves 
D{A){Fig. ^) tied to fiLi/J^n) {Fig. (sji), and the 
unstable mode displays the typical bias depen- 
dent 7®^. The D,A state symmetries implies that 
the mode involves transverse vibrations, see in- 
set of Fig. g In Fig. [4] we show the correla- 
tion between the mter-electrode damping terms 
{^lr)^ and the electron-phonon matrix element 
between the D,A states, {D\M\A). It is seen that 
a larger matrix element leads to higher negative 
m^er-electrode damping terms, while the competi- 
tion with the mtra-electrode damping (and possi- 
bly other damping mechanisms) eventually deter- 
mines if the mode become unstable. We note that 
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the 63meV mode becoming unstable at ~ 0.2V has 
a frequency out side the electrode phonon bands. 
Thus we do not expect additional harmonic damp- 
ing due to these. 




FIG. 4: (Color online) Correlation between the inter- 
electrode damping rate, (incl. only LR terms in 
A, B), and phonon-coupling matrix element between 
the D and A- states. Each dot represents one phonon 
mode. Larger dot has smaller mode frequency, and 
darker one has smaller 7£l + 7rr- Modes displaying 
7^^ < for y < 0.4V are marked with cross line or 
frequency. Insets: The 63meV mode, and bias depen- 
dence of its damping rate. 

Conclusions. — We have discussed how stimu- 
lated amplification of phonons becomes possible 
in molecular rectifiers once the donor level is lifted 
higher than the acceptor by the applied bias. This 
is akin to population inversion leading to 'lasing', 
and may cause instabilities such as switching or 
contact disruption for a certain voltage[4]. On the 
other hand, current-induced cooling may also be 
possible [10 . Both effects draw some analog with 
the current-induced negative damping and cooling 
of a nanomechanical oscillator [30 . 

We thank Prof. Jauho for comments, the 
Lundbeck foundation for financial support (R49- 
A5454), the Danish Center for Scientific Comput- 
ing(DCSC) and Dir. Henriksens Fond for provid- 
ing computer resources. 



* Electronic address: ljtlu@nanotech.dtu.dk| 
[1] A. Aviram and M. A. Ratner, Chem. Phys. Lett. 

29, 277 (1974). 
[2] R. M. Metzger, Acc. Chem. Res. 32, 950 (1999). 
[3] M. Elbing, R. Ochs, M. Koentopp, M. Fischer, 

C. von Hanisch, F. Weigend, F. Evers, H. B. We- 



ber, and M. Mayor, Proc. Nat. Acad. Sci. USA 

102, 8815 (2005). 

[4] I. Diez-Perez, J. Hihath, Y. Lee, L. Yu, 
L. Adamska, M. Kozhushner, I. Oleynik, and 
N. Tao, Nature Chem. 1, 635 (2009). 

[5] N. J. Tao, Nature Nanotech. 1, 173 (2006). 

[6] M. Galperin, M. A. Ratner, A. Nitzan, and 
A. Troisi, Science 319, 1056 (2008). 

[7] J. Hihath, C. Bruot, and N. Tao, ACS Nano 4, 
3823 (2010). 

[8] N. Okabayashi, M. Paulsson, H. Ueba, Y. Konda, 
and T. Komeda, Phys. Rev. Lett. 104, 077801 

(2010) . 

[9] M. Galperin, M. A. Ratner, and A. Nitzan, J. 

Phys.:Cond. Mat. 19, 103201 (2007). 
[10] Z. loffe, T. Shamai, A. Ophir, G. Noy I. Yutsis, 

K. Kfir, O. Cheshnovsky, and Y. Selzer, Nature 

Nanotech. 3, 727 (2008). 
[11] D. R. Ward, D. A. Corley J. M. Tour, and D. Na- 

telson. Nature Nanotech. 6, 33 (2011). 
[12] Z. Huang, F. Chen, R. D'Agosta, P. A. Bennett, 

M. Di Ventra, and N. Tao, Nature Nanotech. 2, 

698 (2007), ISSN 1748-3387. 
[13] R. H. M. Smit, C. Untiedt, and J. M. van Ruiten- 

beek, Nanotechnology 15, S472 (2004). 
[14] G. Schulze et a/., Phys. Rev. Lett. 100, 136801 

(2008). 

[15] M. Tsutsui, Y. Teramae, S. Kurokawa, and 

A. Sakai, Appl. Phys. Lett. 89, 163111 (2006). 
[16] S. J. van der Molen and P. Liljeroth, J. Phys.: 

Cond. Mat. 22, 133001 (2010). 
[17] D. Dundas, E. J. McEniry and T. N. Todorov, 

Nature Nanotech. 4, 99 (2009). 
[18] J.-T. Lii, M. Brandbyge, and P. Hedegard, Nano 

Lett. 10, 1657 (2010). 
[19] B. Persson and M. Persson, Surface Science 97, 

609 (1980). 

[20] M. Head-Gordon and J. Tully J. Chem. Phys. 

103, 10137 (1995). 

[21] T. Frederiksen, M. Paulsson, M. Brandbyge, and 
A.-P. Jauho, Phys. Rev. B 75, 205413 (2007). 

[22] K. Stokbro, J. Taylor, and M. Brandbyge, J. Am. 
Chem. Soc. 125, 3674 (2003). 

[23] We employ the SIESTA,TranSIESTA/Inelastica 
methods [23,21], similar parameters as in Ref. 21, 
and the F-point approximation . 

[24] J. Soler, E. Artacho, J. Gale, A. Garcia, J. Jun- 
quera, P. Ordejon, and D. Sanchez-Portal, J. 
Phys.:Cond. Mat. 14, 2745 (2002). 

[25] L. Ruan, F. Besenbacher, I. Stensgaard, and 
E. Laegsgaard, Phys. Rev. Lett. 70, 4079 (1993). 

[26] T. Nishino, T. Ito, and Y. Umezawa, Proc. Natl. 
Acad. Sci. U. S. A. 102, 5659 (2005). 

[27] G. Schuh, T. Frederiksen, A. Arnau, D. Sanchez- 
Portal, and R. Berndt, Nature Nanotech. 6, 23 

(2011) . 

[28] M. J. Ford, R. C. Hoft, A. M. McDonagh, and 
M. B. Cortie, J. Phys.:Cond. Mat. 20 (2008). 

[29] L. A. Zotti, T. Kirchner, J. C. Cuevas, F. Pauly, 
T. Huhn, E. Scheer, and A. Erbe, Smah 6, 1529 
(2010). 

[30] M. P. Blencowe, J. Imbers, and A. D. Armour, 
New J. Phys. 7, 236 (2005). 



4 



